Enhancement of pair correlation in a one-dimensional hybridization model 
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We propose an integrable model of one-dimensional (ID) interacting electrons coupled with the 
, local orbitals arrayed periodically in the chain. Since the local orbitals are introduced in a way 

that double occupation is forbidden, the model keeps the main feature of the periodic Anderson 
model with an interacting host. For the attractive interaction, it is found that the local orbitals 
enhance the effective mass of the Cooper-pair-like singlets and also the pair correlation in the ground 
state. However, the persistent current is depressed in this case. For the repulsive interaction case, 
the Hamiltonian is non-Hermitian but allows Cooper pair solutions with small momenta, which are 
induced by the hybridization between the extended state and the local orbitals. 



C/3 



0^ 
0^ 



O 
o 



71.10.Pm, 75.20.Hr, 75.30.Mb 



. I. INTRODUCTION 



Metallic compounds containing elements with partially filled d-shells or /-shells, belong to the category of strongly 
correlated electrons. Two typical examples are the high-Tc superconductors and the heavy fermion compounds in 
Q . which the spin fluctuation may play central ruleElo. The normal state properties of heavy fermion compounds are 
Q ' characterized by a large Pauli susceptibility and specific heat as compared to conventional metals. Such phenomena 
are attributed to the large effective mass of the electrons near the Fermi surface. These anomalous are generally 
I ' believed to be due to the formation of resonant states at the Fermi level, which is induced by the admixture of local 
^ , /-orbitals and the conduction electrons, and therefore the systems are usually modeled by the periodic Anderson 
\^ < model or the Kondo lattice model in some limiting cases. One of the majo]j-|mysteries of the heavy fermions is how 
\^ superconductivity could be supported in a system with strong local moment^. It is generally accepted that magnetic 
impurities in BCS superconductors break the time reversal symmetry and are unfavorable to the formation of Cooper 
pairs. Such a pair breaking effect directly causes the reduction of the energy gap of the superconducting state and 
, the transition tcmpcraturetl. However, the situations may be different in some strongly correlated electron systenis, 
where the electrons from the same source could be responsible for both the superconductivity and the magnetismEJ. 
How magnetism and superconductivity reconcile each other is a hot interest in modern condensed matter physics and 
' still remains to be an open problem. 
^ ] Based on the development of the strongly correlated electron systems and low-dimensional systems, niWY efforts 
I ■ have been done in recent years to-|rHjiderstand how the magnetic impurities behave in a ID correlated hostu. Several 
' O integrjable models were proposedB~Q to account for this problem and some novel features were found. In a recent 
^ papercl, Schlottmann studied the attractive Hubbard model with a finite concentration of magnetic impurities. He 
found the impurities generally weaken the binding energy of the singlet pairs and the spin gap could be closed above 
a critical concentration of impurities. p. 
The quantum inverse scattering method (QISM) provides a powerful tool to construct integrable models-Ln IDD. 
In a lattice model, a local operator L„^7.(A) can be defined, which satisfy the following Yang-Baxter relationtj 



X 

I i?r,r'(A, /i)i„,r(A)i„,r'(/-i) = ^n,T'(M)-^ri,r(A)-Rr,r'(-^:M): (1) 

where L„^t-(A) acting on the auxiliary space Vr and the quantum space Vn respectively, A and /i are the spectral 
parameters, i?7-,T'(-^, A*) = ^t,t'(A — /i) is a c- number matrix. Define the transition matrix IV (A) as 

r,(A) =Li,,(A)---Lw,,(A), (2) 

where N is the site number of the lattice. From (1) we can easily show that 

Rr,r'{\ - fl)TriX)Tr,{fl) - T,- (/i)r, (A)i?,,,. (A - /i) . (3) 
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Tracing r and r' in the above equation we have 

[trrTriX),trrTriti)]^0. (4) 
Suppose t(A) = trT-Tr^X) aUows the following expansion 

r(A)=io + | + ^ + ---. (5) 
Since A and fi are arbitrary parameters, from (4) we obtain 

[tm,tn]=0, TO,n=0, 1,2--- (6) 

Choosing one oitn as the Hamiltonian of the system, from (6) we know that all tm are conserved quantities. Therefore, 
we can establish the common eigen states of these quantities. Generally, the generating operators of these eigen states 
are chosen from the off diagonal elements of the matrix Tt-(A). For the impurity models, the impurities are added 
by including some inhomogeneous vertex operators in the transition matrixQ, which satisfy the same Yang-Baxter 
relation to that of the host. 

In this paper, we consider a model of ID interacting electrons coupled with the local orbitals arrayed periodically 
in the chain. Maximumly, only one electron can occupy a single local state. Therefore, the model preserves the main 
feature of a ID periodic Anderson model in the limit U oo. The structure of the present paper is the following: 
In the subsequent section, the model Hamiltonian and its Bethe ansatz solution will be constructed based on QISM. 
In sect. Ill, we discuss the attractive interaction case. It is shown that the effective mass of the Cooper-pair-like 
singlets as well as the pair correlation in the ground state are enhanced by the local orbitals. In sect. IV, we study 
the repulsive interaction case. Sec. V is attributed to the concluding remarks. 



II. THE MODEL AND ITS BETHE ANSATZ 



In this paper, we construct an integrable model which describes interacting conduction electrons in a continuum 
medium hybridizing with local orbitals. The atoms with local orbital are arrayed periodically in the chain and the 
local states are described by the Hubbard operators = |a„ >< /?„|, (a„,/3„ — 0, t,i), with the constraint 

X"-^ +^fj, +^00 — which means double occupation of the same orbital is forbidden. Let us start with the following 
segment transition matrix 

dxTn+i{X\x,na) =: L{X,x)Tn+i{X\x,na) :, (7) 

where x € [na, {n+ l)a); :: denotes the normal order of the fermions; a is the space between two nearest local orbitals; 
the L operator is defined as 

/ fA iVffCt(x) 

with cj.(a;) {ccr{x)) the creation (annihilation) operator of the conduction electrons. Here we take the boundary 
condition of Tn+i{X\x,na) as 

/a'(A)-&'(A)Xf^ -&'(A)Xp^ b'{X)X^^ 
T„+i(A|na,na) = L„(A)= -b'{X)X^^' a'(A) - 6'(A)X«^ fe'(A)Xf\ 

V 6'(A)Xo"; &'(A)Xo"_x ' a'(A) + 6'(A)X«o 

where a'(A) = a(A -I- i.g/2), fe'(A) = 6(A + ig/2) and 

«(A) = biX) = (8) 

X — tg X — ig 



We remark that with an unit, boundary condition, the transition matrix is just that of the (5-potential Fermi gas 
modeltJ introducedpby Yangt2l. The non-unit boundary condition is much similar to the inhomogeneous L-operator 



rqmar 

„„^_eeii , 

in the lattice modelEl. It is easy to show that the following Yang-Baxter relations hold 



2 



R{X - fi)L{\, x) (g)s L{fi, x) = L{fi, x) (g)s L{X, x)R{X - fj.), (9) 

R{X - fi)L„{X) (g)s Lnifi) = i„(^) (g), Ln{X)R{X - n), (10) 

where 

R{X) =a{X)P + b{X), (11) 

and P is the FFB graded exchange operator acting on the direct product of the auxihary spaces, -P^^^^ = 
<5oi&2^a26i (— l)*^'^*^'^ ; ^1,1 = IjEo = (for convpiience, we put the subscripts t,i and as 1,2 and 3 respectively); (g)^ 



denotes the direct product with FFB gradingt3llil 

{F ®, G)fa = FabG,d{-ir'^''+''^. (12) 

From (9) and (10) we can easily derive 

R{X - fi)Tn{X) Tnifi) = T„(/i) ®, T„(A)i?(A - m), (13) 

with T„(A) = T{X\na,na — a). Furthermore, the global transition matrix 

T{X) = Tn{X)Tn-i{X)---Ti{X) (14) 

satisfies the same Yang-Baxter equation of the segment ones T'„(A), 

RiX - Ai)r(A) ®s Till) = Tip) ®, TiX)RiX - /^), (15) 

where N is the number of atoms with local orbital. Introduce the notation 

/An (A) Ai2(A) Bi(A) 
T(A) = A2i(A) A22(A) B^iX) 
\ Ci(A) C2(A) DiX) 

From (15) we have the following commutation relations 

AabiX)C,iti) = (-l)^»'.!l^Llilltcp(M)A,,(A) + ^[^^ ^\ aiX)AM, (16) 

a(A — /i) a(A — /i) 

^^^^^^^''^ ^ ^ii^^f'^^^^^^^ ~ ^(^f'^^^^^^"^' ^^^^ 
a,(Ai)a,(A2) =r(Ai - A2)^j:ja,(A2)a,(Ai), (18) 
[T(A),T(/i)] =0, (19) 

where 

t(A) = strT(A) = -Aii(A) - A22(A) + DiX) (20) 

and 

ri^Tcd = ~biX)dabScd - a(A)5ad'5bc. (21) 
It can be shown that r(A) satisfies the following Yang-Baxter relation 

'■(A - ^^)ZZ<^)T2ri^.tfH = rii,):iZriX)ZilriX - /.)-^^. (22) 

From (19) we can see r(A) can be treated as a generator of an infinite number of conserved quantities. Choose the 
vacuum state as 

c,(a;)|0>=Xo"_,|0>=0. (23) 
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We have 

/e'^^a^{X + q) 

T(A)|0>= e'^^a^{X + i§) | |0 > . 

V Ci(A) C2{X) e-'i^ 

Therefore, Ca(A) can be treated as the creation operators of the eigenstates of r(A): 

\h,---,kn\F >= (fci)C„, (fc2) • • • C„„ (fc„)|0 > , (24) 

where the indices Uj run over the vahics 1,2 and i^""'""! is a function of the spectral parameters kj. Prom the 
commutation relations (16) and (17) we have 

D{X)\ki,---,kn\F>=e-'i'^r\— -\ku---,kn\F> 

n n 

+ E(^OaV.-.-t f^;-. W n > ' (25) 

n n 

[An (A) + A22{X)]\ku- ■ ■ , kn\F >= -e'^^a^(A + i|) J] _, . H ^^-^ > 

x[r(^HA)]^;;;t^"''---"^ +E(A;)^rtc'6,(A)nc,,(fc,^^ > (26) 

1=1 m 



where 



and 



t(i)(A) = str[T^^\X)] = str[L<<^\X - fc„) • • • ^(^^(A -k,)], (27) 
LW(A) = 6(A)pW+a(A) (28) 

which fulfill the Yang-Baxter relation 

r(A - H)L^\X) ®, X«(m) = 4'^(/.) ®. 4'^(A)r(A - /.), (29) 
with (-P^^^)cd = —^ad^bc the 4x4 permutation matrix; 

a(A - fc;) a(A;; - fc,) 2' 
x(-l)'+i4^)(fc; - fc„)^lt-L«i(fc, - kr^-iizltz: ■ ■ ■ Li+i{ki - k+i)tZa,^^, (30) 



^ ' ' ^WV...a,^ [ a{ki-Xy lLa{ki-ki)' ^ ^ 



where 



S{ki)t-t = - kt'-\Z-. ■ ■ ■ <ki - ki^. (32) 

To obtain the Bethe ansatz equations and the eigen values of r(A), we solve first the following eigen value problem 

T^^\X)F = e{X)F. (33) 

Prom (29) we know 
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r(A - f^)Ti'\X) ®, t(i)(m) = Ti^)(M) ®. T(i)(A)r(A - „). (34) 

Introduce the notation 

(A) - ^c^D^x) DW{X) ) ■ 

(34) gives 



Acting T^i^A) = -A(^\X) - D^^\X) on (38), we have 



M n , s M 



and the cancellation of the unwanted terms gives the Bethe ansatz equation 

n M , ■. 



(35) 



A(i)(A)C(i)(M) = ^^^C(i)(M)A(i)(A) + ^|^C7(i)(A)^(i)(M), (36) 

C(i)(A)C(i)(m) = C(i)(/x)C(i)(A), [r(i)(A)(A),r(i)(M)] = 0. (37) 
Define the pseudo vacuum |0 >(i) as S(i)(A)|0 >(i)= 0. The eigenstates of t(i)(A) can be written as 

M 

\i^i,---,m>=l[c^'\i^a)\o>^^^ . (38) 



(40) 



To ensure (24) to be an eigenstate of r(A), the unwanted terms in (25) and (26) must cancel, i.e., 

+ (AO^lV-tli^""-"^ = 0. (41) 

This gives 

M 

a^(fc,+z|)e*'='^ = I]a(fci-Ma)- (42) 



(For the detailed derivation of (30), (31) and (42), we refer the readers to the appendix B of Ref.[12], while the 
algebraic structure of the present model is almost the same to that of theirs.) The eigenvalue of t(A) reads 

.(A, W}) = a-(A + .|)e«^- ^^fW + e^^^ ^(^' (43) 

where Ng = nis the electron number and M is the number of electron with down spins. Putting ^,a = Xa — ig/2, the 
Bethe ansatz equations (40) and (42) are reduced to 



N 



M 



^ikiNa I "'J ' 2 

e 



kj + -kg \ -TT kj Xc, + ig 



n Z_TLZ ' (44) 



-TT Aq, - Kj - ^g ^ _ TT Ag - A/3 - ig 
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and the eigenvalue of r(A) can be written clearly as 

We 



,y{X, {A4) = e-^^^ + 

1 = 1 3 



-6 2 



^ ^ A + i| ^ I- -p-|- A — Aq — i| Yl l^ ^9 \ , TT n I 



Now we turn to the construction of the model Hamiltonian. For A — > oo in the upper half complex plane, we have 
the following asymptotic expansion 

ln[r(A)e^*^] = 1 + ^g{^ + ^ + ^ + . . .}. (47) 

From the commutation relation (19) we know 

[a„,C„]=0, m,7i= 1,2,3, (48) 

In principle, we have more freedoms to choose a Hamiltonian from the conserved quantities {C„}. In this paper, we 
define the Hamiltonian as 

H = C3+t9C2-^9^C,. (49) 
For imit boundary condition Tn+i{na, no) = 1, (49) can be derived via Neumann expansion aJl3 

^o-E/ + / c\{x)cl{x)c^{x)c^{x)dx, (50) 

which is nothing but the Hamiltonian of a ID electron gas with (5-potential interactions. With the local orbitals, (49) 
reads 

H = Ho + Hj, (51) 

where Hj is very complicated, which contains a hybridization term J2n k^n^kcl,{k)X^ ^ + h.c, a correlation term 
X]„o-(t' Ua,a'cX{na)c„{na)X^, and other irrelevant terms which often appear in the integrable impurity modelsB. 
Comparing (47) and (46), we easily obtain the spectrum of the Hamiltonian is given by 

E = Y.kl (52) 



III. ATTRACTIVE INTERACTION CASE 



We discuss first the attractive interaction, i.e., 5 < case. Without the local orbitals, the ground state is a Fermi 
sea filled by the Cooper-pair-like bound pairs. Now by studying (44) and (45), we can show how the local levels behave 
in the ground state. Carefully checking the Bethe ansatz equations we find the the Cooper pair states described by 

fc± = A„ ± ^5, (53) 

are still possible solutions in the thermodynamic limit L ~f 00, despite the existence of the local orbitals, where Aq 
are real. To study the stability of these pair states, we consider a reference state, i.e., all the A^e electrons form Ne/2 
pairs. In this case, the Bethe ansatz equations are reduced to 



Aa+^lffiy A„-A;3+*|g|' ^ > 
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and the energy of this state is given by 



E=Y.^>^1-^- (55) 



C( = l 



Taking the logarithm of (54) we obtain 



1 ^^'^ 

2K + -e{K)^—^+y^e{K-\p), (56) 



where 9(x) = 2ia,'n~'^{x /\g\) and are integers or half integers depending on the parity of iV — Ne/2. Notice that 
each /q, corresponds a pair state and they must be different from each other due to the exclusion principle. The 
minimum state (with lowest energy) is thus described by a sequence of {la} = {—[Ne/2 — l)/2, • • • , {Ne/2 — l)/2}. 
In the thermodynamic limit Ne/ L = n ^ finite, the distribution of Aq, can be described by a density fmiction 

p{Xa) = lim — ^— (57) 



which satisfies 



TT a 

where the cutoff Af is given by 



p{X) = 1 + i/(A) - /(A - A'))p(A')A' (58) 

J-Xf 



n 

p{X)d\ = -, (59) 



and /(A) ~ |.g|/7r(A^ + g^). Without the local orbitals, the density distribution of A in the ground state takes the form 



1 '■^ 



TT 



Po(A)---/ /(A-A')po(A')dA' (60) 



77 

Po(A)dA = -. (61) 

AO, 2 

Comparing (58) and (60) we can readily read off Xp < X'^. That means the effective Fermi energy is reduced by the 
hybridization relative to that of the homogeneous system, a typical heavy fermion behaviom. Now we consider the 
density of states at the Fermi surface. The energy density of the minimum state (relative to the chemical potential) 
can be written as 

E/L - (2A2 - ^ - 2fi)p{X)dX, (62) 

where p is the chemical potential. Substituting (58) into (62) we readily get 

E/L^ r [- + -f{XMX)dX, (63) 

where e(A) is the dressed energyEl of the Cooper pairs 

e(A) = 2A2 - ^ - 2/i - /(A - A')e(A')dA'. (64) 

Consider a particle-hole excitation relative to the minimum state. The density of the hole and that of the excited 
particle can be expressed as 

p,(A) = i,5(A-A,0, Pp(A) = l<5(A-Ap), (65) 
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where Xh and Xp are the centers of the hole and the particle, respectively. In addition, such an excitation induces the 
back flow of the Fermi sea, i.e., Sp(X). From the BAE we have 

p(A) + SpiX) = - + -fix) + p„(A) - pp{X) 
TT a 

fix - x'Mx') + dpix') + ppix') - ph{x')]dy. (66) 

With (58) we have 

6p{X) = ph{X) - ppiX) + ^[/(A - Xp) - fix - Xi,)] - fix - X')SpiX')dX'. (67) 
The excitation energy reads 

6E = L [ 2X^SpiX)dX-2Xl + 2Xl. (68) 

Substituting (67) into (68) we readily obtain that 

SE = eiXp)-eiXh). (69) 

Therefore, e(A) can be treated as the quasi-particle energy of the elementary excitations. The excitation of breaking 
a Cooper pair can be treated in a similar process. Such an excitation can be described by a A-hole in the Fermi sea 
and two real A;- modes ki , k2 above the Fermi level. In this case, the excitation energy is 

5E = -eiXh)^2p + kj + kl + ^. (70) 

Notice the dressed energy has the properties e(±AF) = 0, e(A) < for |A| < A_f and e(A) > for |A| > A_f. Concerning 
the excitation near the Fermi surface, i.e., Xh — > Xp, fci,fc| — > p, from (70) we readily obtain that there is a finite 
gap to break a Cooper-pair at the Fermi surface. The energy gap seems not changed via the local orbitals. 
The minimum state we introduced is thus the absolute ground state of the system for a given N^. From the BAE 
(56) we can see the "quasiparticle" (Cooper pair) momenta can be defined as p(Aq) = 2TrIa/L, then in our case 
p'iX) = 2-KpiX). The density of states at the Fermi surface is thus 

IdpiX) _ 2piXF) _ 1 

- ^MX)^'='^ - Iv^ = ^' 

where v is the sound velocit Since piXp) > PoiXp) and e'iXp) is an increasing function of A^ as can be shown 

in (64), we deduce that the density of states is enlarged by the local orbitals. It is not very strange because the Fermi 
sphere is compressed. 

At present, we can see that the local orbitals can not destroy the Cooper pair state completely. However, it is still 
not clear whether the local orbitals weaken the pair correlation or enhance it. To answer this question, let us consider 
the stiffness constant K which measures the non-universal exponents of apvariety of correlation functions in IDtj. For 
the integrable systems, K = Z^(A) and the dressed charge function Z(A)ll3 in our case satisfies 

Z(A) = 1- f " fix- X')ZiX')dX'. (72) 

Notice the local orbitals do not change the form of Z(A) but the value of the cutoff Xp. Easily we can show 
dZiX) / dXp < 0. That means ZiX) is a monotonically decreasing function of Xp. The stiffness constant K is therefore 
also a monotonically decreasing function of A. As stated in earlier publicationsli3, the gapless ID quantum system 
is conformally invariant at zero temperature and the non-uniYersal exponents of the correlation functions can be 
derived from the finite size corrections of the energy spectrural3'E3: 6E = 2iTvXy/ L, where Xy is the scaling dimension 
(one half the the critical exponent) of the relevant operator. In our case, the spin excitations have a finite gap 
while the charge ones are gapless. Therefore, the charge sector is conformally invariant at zero temperature and the 
asymptoticlong-distance superconducting correlators can be derived from the finite size correction of the ground state 
energyOotj. Notice the pair operator cj(a::)c|(a;) induces a pair number change by one. The energy change induced 
by this operator can be calculated by following the standard method introduced in Ref.[13] as 
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SE^Ei,^^-Ei,-2, = i^^=J^^. ,73) 

where is the ground-state energy with n Cooper pairs. Hence the pair correlator reads 

< c^{x)ci{x),c\{0)c\{0) >^ X'' (74) 

Since 9 < Oq (where 6o is the corresponding exponent of the homogeneous system) , we conclude that the local orbitals 
enhance the superconducting correlations. To see it clearly, let us consider the tunneling of the Cooper pairs through 
an impurity. The leading tunneling current through the impurity is 

J{x,y) i[cl{x)c\{x)c'f{y)ci{y) - h.c], x 0~ , y ^ , (75) 



where we have put the impurity at the origin. From the boundary conformal field theoryllZl we have the time correlator 
of J 

[j{t),.m]^t-^'^. (76) 

The tunneling conductance can be derived from the Kubo's formulaO as 

(t(T) - T^^"'"^ (77) 

That means the local orbitals enhance the low-temperature tunneling conductance, which provides another evidence 
of the enhancement of the superconducting correlation at low temperatures. The Drude weight D, which measures 
the persistent current in lDt3, can also be derived exactly for the present model. It reads D = Kv. From a similar 
discussion, we can deduce D is reduced by the local orbitals. That means the local orbitals depress the persistent 
current. 

IV. REPULSIVE INTERACTION CASE 

Now we turn to the repulsive case {g > 0). In this case, the Hamiltonian defined in (49) is non-Hermitian. We 
remark the study on non-Hermitian systems has drawn considerable attention recently for ptiie applications in a 
variety of physical situations sucK-as delocalization in the disordered systemala, quantum chaosESI, population-biology 
in random media with convectiorE^ and metal-insulator transitions derived by|-an imaginary vector potentialE3. The 
spectrum of a non-Hermitian Hamiltonian generally falls in the complex planec3. In our case, from the Bethe ansatz 
equation (44) we can see the solutions of kj are classified into two types: (i)|0(fcj)| = 1; (ii)|0(fcj)| 1, where 

^(fc)=e^'=-:^±i^. (78) 

For the Hermitian Hamiltonian, case (i) represents the real modes and case (ii) denotes the string solutions (Cooper 
pairs in the present model). Both of these two types of modes give out real eigenvalues of the Hamiltonian. When 
g > in our case, the Hamiltonian contains a non-Hermitian hitbridization term and a Hermitian term which is 
nothing but the Hamiltonian of the repulsive (5-potential Fermi gaslliJ. Even for |0(fc)| = 1, the Bethe ansatz equations 
have complex solutions. A typical solution is the imaginary mode k = in with na = ln[|«; -I- g/2\/\K — g/2\ in the 
thermodynamic limit. An interesting feature is that the Bethe ansatz equations allow Cooper pair solutions. For 
|(/)(A:j)| > 1, the left hand side of (44) is divergent in the thermodynamic limit while for |(?I)(fcj)| < 1 it tends to zero. 
Therefore, the pair solutions (53) are allowed. However, there is a constraint for the real part of the pair solutions. 
The condition \4>{k)\ ^ 1 hints 

e-*"|l + ^|>l, (79) 
A 

which gives the critical value of the real part of the pair solutions 

Ac = , ^ ■ (80) 
Vef" - 1 

The Cooper pair solution can exist only in the region —Ac < A < Ac- For the complex solutions in case(i), |/mfc| > \g\/2 
when \Rek\ < Ac. Hence the ground state consists of these modes and there is no Cooper pair in the ground state. 
Even so, we can see the local orbitals do cause some fluctuation towards the formation of Cooper pairs, since without 
these orbitals, the system does not have any bound state of electrons. 
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V. CONCLUSION 



In conclusion, we propose an integrable hybridization model for a ID correlated electron system. We remark 
though the local orbitals are introduced periodically, the problem is still at the level of single impurity since there 
is no correlation among the local states. In fact, the Bethe ansatz equations (44) and (45) do not depend on the 
distribution of the local orbitals in the real space. It remains still a shortcoming of the integrable models with many 
impurities since the impurities introduced in such a way are completely transparent to the host electrons and only 
the forward scatterings are included. Even so, a finite concentration of impurities enhances the superconducting 
correlation in our model. It seems that our result contradicts to that of Rcf. [8]. The difference comes from (i) in 
our model, the binding energy is independent of the momenta of the pairs while in the attractive Hubbard model, it 
depends on the pair momenta (A in [8]); (ii) there is no Eimp in the formal expression of the eigenenergy (52), though 
the local orbitals change the distribution of {kj}; (iii) in our model, the density of electrons is not changed by the local 
orbitals. This corresponds to the hybridization configurations of /° ^ for the local orbitals with a configuration 
/"^ in the atomic limit. We remark impurities with higher spin and finite momenta can also be introduced in our 
model with a similar procedure. In this case, the spin momenta are split into two classes. One associated with the 
pair states (53) and the other represents the dynamics of the rest spin degrees of freedom of the local orbitals. The 
latter excitations are always gapless as in the spin chains and the spin excitations breaking a pair are still gapful. 

One of the authors (YW) acknowledges the financial supports of Alexander von Humboldt-Stiftung and China 
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